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The troposphere is the part of the atmosphere where people live. The troposphere goes 
up to about 12 km above the earth over places at middle latitudes like Washington, D.C. 
The next layer of air up in the atmosphere is called the stratosphere. The tropopause 
separates the troposphere and stratosphere. Most of the ozone in the atmosphere is in the 
stratosphere where it protects people from the harmful rays in sunlight. Near the earth’s 
surface, ozone is a pollutant. In general, air travels from the troposphere to the 
stratosphere through upward motion in the tropics. Most of the air comes back down to 
the troposphere at middle latitudes. The so-called “stratospheric residual circulation” 
moves the air in the stratosphere from the tropics to the middle latitudes. In this way, the 
amount of air that moves out of the stratosphere into the troposphere at middle latitudes 
depends on the strength of the residual circulation. This overturning of the atmosphere is 
important as it brings chemicals produced by the activities of people into the stratosphere. 
Most notable are chlorofluorocarbons (CFCs) that destroy stratospheric ozone. The 
transport of air out of the stratosphere removes ozone-destroying chemicals, but also 
brings ozone into the troposphere. 

Previous studies using atmospheric models have shown that the residual circulation and 
the amount of air that crosses the tropopause change because of increased chemicals 
produced by people in the atmosphere. The strength of the residual circulation will 
contribute to how much of a chemical can be found in the stratosphere at any time. 
Increasing the speed of transport from the stratosphere to the troposphere would also 
increase how fast ozone-destroying chemicals are removed from the stratosphere. Thus, 
knowing the strength of the residual circulation is needed to predict when the amount of 
ozone in the stratosphere will return to the greater levels that existed in the past. 

This work shows a method of detecting a change in the residual circulation from 
observations of hydrogen chloride (HC1). HC1 in the atmosphere is primarily produced 
by chemical reactions in the stratosphere. For several years, instruments on the ground 
have made measurements of how much HC1 exists in the atmosphere directly above 
them. This type of measurement is called the “column HC1”. The middle latitude 
column HC1 tends to be greatest in the late winter and least in the late summer. The 
amount that the column HC1 changes within a year tells us how strong the residual 
circulation is. Many measurements are needed each year to find out how much the 
column HC1 changes. So far, not enough years have enough measurements to show if the 
residual circulation is changing. This work shows that many instruments working 
together for a few more years could determine if the residual circulation is changing. 
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Abstract 

Ground-based measurements show that the HC1 column at middle latitudes varies 
seasonally with a late spring maximum and a minimum in late summer. HALOE 
measurements of HC1 profiles show that this seasonal cycle in the column is controlled 
by the seasonal variations of the mass and composition of the lowermost stratosphere. 
The seasonal maxima are produced by wintertime downward transport of stratospheric 
air. We hypothesize that the amplitude of the seasonal cycle is quantitatively related to 
the strength of the stratospheric residual circulation and stratosphere-troposphere 
exchange (STE). We investigate this relationship using a multi-decadal simulation from 
the GSFC Chemistry and Transport Model driven by output from a general circulation 
model in which an STE trend exists. We also examine the ground-based record of HC1 
measurements at Jungfraujoch for evidence of a trend in the HC1 seasonal amplitude that 
would be a signature of a long-term change in the stratospheric circulation. 
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1. Introduction 


General circulation models have been shown to exhibit a change in the residual 
circulation as a result of climate-induced change in atmospheric composition [e.g., Rind 
et ah, 1998; Sigmond et al., 2004; Eichelberger and Hartmann, 2005]. If such a trend is 
real, it would contribute to observed trends in stratospheric trace gases including ozone. 
Butchart and Scaife [2001] have demonstrated that increasing the rate of mass exchange 
between the stratosphere and troposphere accelerates the removal of chlorofluorocarbons 
(CFCs) from the stratosphere and would alter ozone recovery. The present study 
introduces a novel method of detecting a change in the residual circulation. 

Both the change in the mass of the lowermost stratosphere and the change in the HC1 
mixing ratio dominate the seasonal variation observed in the column. The mass of the 
extratropical lowermost stratosphere is a maximum in the winter when the tropopause is 
lowest and a minimum during summer when the tropopause is highest [Schoeberl, 2004]. 
The composition of the middle latitude lowermost stratosphere is dominated by 
downward transport of stratospheric air from the region above 380K in the winter. In the 
summer, sub-tropical tropospheric air is quasi-isentropically transported into the middle 
latitude lowermost stratosphere (the region between 380K and the tropopause), resulting 
in a composition that is more tropospheric in character [e.g., Ray et ah, 1999; Pan et ah, 
2004]. Middle latitude ground based measurements of the HC1 column vary from a 
maximum in late winter to a minimum in late summer. The stratospheric column above 
100 hPa from profiles measured by the Halogen Occultation Experiment (HALOE) on 
the Upper Atmosphere Research Satellite are nearly constant throughout the year, thus 
the seasonal variability of the lowermost stratospheric mass and HC1 are manifest in the 
column measurements. 

Based on the discussion above it seems likely that the amplitude of the seasonal cycle of 
mid-latitude column HC1 can be used as a marker of the residual circulation, The annual 
maximum in the HC1 column is sensitive to a change in the stratospheric residual 
circulation through the downward flux of HC1 into the lowermost stratosphere. A 
stronger (weaker) residual circulation will produce a larger (smaller) late winter 
maximum in the HC1 column. The seasonal minimum of the column will remain 
relatively unchanged since it is caused by the flux of low HC1 tropospheric air into the 
lowermost stratosphere. Thus, greater seasonal amplitude of the HC1 column will exist 
under the condition of a more vigorous residual circulation. Here we use a multi-decadal 
simulation with a trend in the mass stratosphere-troposphere exchange (STE) to quantify 
the relationship between STE and the mid-latitude column HC1 annual amplitude. We 
also examine the existing ground-based observation record for such a signature of an 
increasing residual circulation. 

2. Simulation and Data 

(a) Chemistry and transport model and simulation 
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We performed a simulation for 1973-2023 using the GSFC Chemistry and Transport 
Model (CTM) [Stolarski et al., 2006]. Boundary conditions of source gases, including 
observed and projected CFCs and other sources of chlorine, are taken from scenario A2 
[WMO, 2003]; Stolarski et al. [2006] provide details of the simulation. Meteorological 
fields from the GEOS-4 general circulation model (GCM) [Bloom et al., 2005] are used 
to drive the CTM. The GCM uses Hadley Centre monthly mean sea surface temperatures 
for the years 1949 to 1998. Douglass etal. [2003] use comparisons of observations of 
long-lived tracers with fields from a CTM simulation driven by an earlier version of 
GEOS-4 to show that the large-scale stratospheric transport of the GCM is realistic. 
Schoeberl et al. [2003] show that the age-of-air and residual circulation produced by 
these GCM fields are reasonable. Olsen et al. [2004] show that the flux of mass and 
ozone from the stratosphere to the troposphere compares well with estimates based on 
observations. 

(b) HALOE and FTIR observations of HCl 

The Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research 
Satellite (UARS) measured die composition of the stratosphere and mesosphere by solar 
occultation limb sounding from 1991-2005 [Russell et al., 1993]. Thirty vertical profiles 
of seven chemical species, including HCl, were made daily. In 1997, power sharing 
necessitated by battery problems reduced HALOE coverage, 


Regular ground based measurements of column HCl have been made by Fourier, 
Transform Infrared (FTIR) spectrometers of the Network for the Detection of 
Stratospheric Change (NDSC). The longest on-going record exists at Jungfraujoch 
extending back to 1989. The random and systematic errors of each measurement are 
±3% and ±5%, respectively [Rinsland et al., 2003]. 


3. Results 

The simulated stratospheric HCl reproduces the main features of HALOE observations 
(Figure 1). Similar to the observations, the seasonal cycle in the mid- to upper- 
stratosphere is relatively flat and almost all of the variability in the column of HCl comes 
from the lowermost stratosphere. The simulation is biased high with respect to the 
observations but is within the calibration error estimates of HALOE [Russell et al., 1996]. 
Any mean bias of the simulated HCl compared with observations is irrelevant to this 
study and is removed from the comparison in Figure 1. The key point of this comparison 
is that the seasonal cycle in the lowermost stratosphere is reproduced well by the 
simulation and the seasonal amplitudes are similar. 

Figure 2 compares the monthly mean Fourier Transform Infrared Spectrometer (FTIR) 
HCl column measured at Jungfraujoch (46.55° N) with the simulated column at 50° N. 
We choose this comparison with a latitudinal offset (one model grid point) so that the 
mean magnitudes are similar. This offset does not alter our conclusions presented in this 
work. Overall, the interannual trend and seasonal amplitude of the simulated HCl 
column compare well with the available FTIR data. The interannual trend reflects the 
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specified boundary conditions of chlorine species. The CTM simulation includes the 
aerosol loading of two major volcanic eruptions, El Chichon in 1982 and Pinatubo in 
1991. The eruption of Pinatubo in 1991 results in a substantial decrease of the seasonal 
amplitude of the zonal mean simulated HC1 column. In comparison, significantly low 
values of column HC1 are observed over Jungfraujoch during winter 1991/1992. 

We have calculated the annual extratropical STE of mass in the GCM simulation using 
the mass balance diagnostic of Olsen et al [2004]. The STE exhibits an increasing trend 
of 3% per decade in the Northern Hemisphere (Figure 3) and 2% per decade in the 
Southern Hemisphere (not shown), This trend in the mass STE and the forcing will be 
discussed in a future work. The focus of the present study is the relationship of the 
seasonal amplitude of the HC1 column to STE. 

We determine the zonal mean column HCl seasonal cycle amplitude by fitting a linear 
function of an annual sine and cosine. The amplitude is divided by the annual mean and 
expressed as a percent to remove the trend of atmospheric chlorine. This normalization 
also removes the impact of the volcanic aerosol loading. 

The time series of the zonal mean column HCl seasonal amplitude at 50° N from the 
simulation is shown in Figure 4a. Over the entire time series, the HCl amplitude 
increases at 9.7% per decade. We are currently preparing a paper that shows that a 
change in the STE is coupled with a change in the stratospheric residual circulation. 
Figure 4b shows that the simulated annual HCl amplitude and STE are significantly 
correlated at greater than 99% confidence, consistent with the trend in the residual 
circulation as the cause of the increasing trend in the HCl seasonal amplitude. The 
residual circulation driven flux of HCl from the region above 380K to the lowermost 
stratosphere increases, and this results in greater column HCl seasonal amplitudes. 

To separate the change in residual circulation from the changes in chlorine loading, a 
second simulation for 1979-2001 was run with source gas boundary conditions fixed to 
1979 levels. All other boundary conditions and driving meteorological fields were 
identical to the longer simulation. The column HCl seasonal amplitude increases at a rate 
of about 13% per decade and is also significantly correlated with STE in the fixed 
chlorine simulation (not shown). Over the same time period, the varying chlorine 
simulation exhibits a 20% per decade trend. Thus, the strength of the residual circulation 
is a primary factor in determining the amplitude of the column seasonal cycle. However, 
the trend in atmospheric chlorine still contributes to the HCl column amplitude trend with 
our normalization by the annual mean. It remains advantageous to use this normalization 
instead of normalizing by atmospheric chlorine since aerosol influences are eliminated 
and it can be easily applied to ground-based data records. Atmospheric chlorine is 
projected to continue to slowly decrease at present day rates over the next few decades, 
so the significant correlation of the HCl amplitude and STE will persist. 

The HCl seasonal amplitude has considerable inter-annual variability. We can estimate 
the number of years of observations that are needed to isolate a trend in the residual 
circulation given the same trend and standard deviation of the HCl amplitude as found in 
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the simulation. Our results from the time series of the HC1 amplitude at 50° N (Fig. 4a) 
indicate that 1 1 and 19 years are needed to isolate the signal with 95% and 99% 
confidence, respectively. If the actual trend of the amplitude is less or if the variance is 
greater, then the number of years needed will increase, and vice versa. The results using 
time series at different latitudes within 38° N to 50° N do not vary greatly. At these 
latitudes the calculated number of years needed range from 10 to 13 years at the 95% 
confidence level. 

Figure 5 shows the seasonal amplitudes calculated from the FTIR measurements at 
Jungfraujoch. We use a Monte Carlo bootstrap method following Hood et al. [1993] to 
estimate the error bars for the fit of the annual sine and cosine. There are many 
significant gaps in the early years of the data record as can be seen in Figure 2. Fewer 
than seventy days have measurements in each of the years 1989-1994 and many months 
lack any measurements. These years are significantly under-sampled, resulting in very 
large estimated errors in the seasonal amplitude. Beginning in 1997, at least 100 
measurements per year are recorded and all months are represented. A positive trend 
exists after the mid 1990’s, but it is not significantly different from zero given the 
variance and length of the record. Thus, no evidence of an increasing residual circulation 
can be determined from the present data. 

In the simulation, we used zonal mean HC1 columns at a single latitude to compute the 
seasonal amplitude. The variability of the column HC1 is much greater at a single 
location than the zonal mean. However, the number of sites of ground-based 
measurements of column HC1 from the NDSC is few. We have found that the mean 
measurement from several middle latitude sites reduces the single station geographical 
and meteorological variability. We have simulated NDSC measurements by sampling the 
model at grid points nearest to Northern Hemisphere middle latitude NDSC sites that 
have any multi-year record of FTIR measurements. The variability of the average of five 
locations between 38° and 47° N (Jungfraujoch at 8° E, Zugspitze at 1 1° E, Moshiri at 
142° E, Rikubetsu at 144° E, and Mt. Barcroft at 1 18° W) is similar to the zonal mean. 
Using the mean of the values sampled at these five sites we find that 1 1 years are needed 
to detect a signal at 95% confidence, the same length of time estimated for the zonal 
mean. This suggests that coordination of simultaneous mid-latitude NDSC FTIR 
measurements to reduce longitudinal variability could identify a trend in a reasonable 
time frame. 

4. Discussion and Conclusions 

Of the trace gases we investigated, ozone, methane and CFC’s, HC1 is most sensitive to 
the residual circulation and STE since it has a strong gradient in the lowermost 
stratosphere and the seasonal cycle above 380K is flat. The seasonal cycle of ozone in 
the lowermost stratosphere is a convolution of transport driven changes and the seasonal 
variation in ozone above 380 K [e.g., Olsen et al., 2004]. The seasonal cycle of other 
long-lived gases such as methane and fluorocarbons are not as sensitive to the residual 
circulation since their vertical gradient in the lowermost stratosphere is relatively small. 
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In summary, we have identified a method of detecting a change in the residual circulation 
from observations of HC1. The seasonal amplitude of the middle latitude HC1 column 
provides a marker of the magnitude of the residual circulation once the interannual trend 
of atmospheric chlorine is removed. Currently, a well-sampled record of ground-based 
observations is not long enough to confidently isolate a signal of increased circulation. 
Our results indicate that such a signal could be detected in a reasonable time with 
coordination by mid-latitude NDSC sites. 
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Figure 1 . Seasonal cycle of partial pressure of HC1 from (a) HALOE at 45° - 55° N and 
(b) CTM at 46° - 54° N. The annual mean at each pressure surface has been removed to 
emphasize the seasonal cycle. 



Figure 2. Time series of HC1 column. Bold line is the model zonal mean at 50° N. Thin 
red line is the monthly mean FTIR data from Jungfraujoch (46.55° N). 
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Figure 3. Annual extratropical STE for the Northern Hemisphere from a fifty-year 
integration of the GEOS-4 GCM. Negative values indicate net transport from die 
stratosphere to the troposphere. Values are smoothed using a five-year running mean. 
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Figure 4. (a) Time series of zonal mean HC1 seasonal amplitude from CTM at 50° N. 
(b) Correlation of seasonal amplitude with annual extratropical STE. 
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Figure 5. Time series of the seasonal amplitude of column HC1 measured by the FTIR at 
Jungfraujoch. 
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